Circuit Theory for Full Counting Statistics in Multi- Terminal Circuits 
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We propose a theory that treats the current, noise, and, generally, the full current statistics of 
electron transfer in a mesoscopic system in a unified, simple and efficient way. The theory appears 
to be a circuit theory of 2 x 2 matrices associated with Keldysh Green functions. We illustrate the 
theory by considering the big fluctuations of currents in various three-terminal circuits. 
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The field of quantum noise in mesoscopic systems has 
been exploded during the last decade, most achievements 
being summarized in a recent review article. Q Measure- 
ment of fractional charge in Quantum Hall regime [Q, 
noise measurements in atomic-size junctions || and su- 
perconductors Q] are milestones of the field and demon- 
strate the importance of quantum noise as a unique tool 
to study electron correlations and entanglements of dif- 
ferent kinds. A very important step has been made in |j| 
where an elegant theory of full counting statistics (FCS) 
has been presented. This theory encompasses not only 
noise, but all higher momenta of the charge transfer. 

Starting from pioneering work of BiittikerjfJ, a special 
attention has been paid to noise and statistics of elec- 
tron transfer in multi-terminal circuits. The correlations 
of currents flowing to different terminals reveal Fermi 
statistics of electrons. These cross-correlations have been 
recently observed. Although the noise correlations for 
several relevant layouts have been understood 0], the 
evaluation of FCS still encountered difficulties. For in- 
stance, an attempt to build up FCS with "minimal corre- 
lation approach" || has lead to contradictions . || This is 
unfortunate, since higher-order current correlations sup- 
ply information about higher-order electron correlations 
and multi-particle interference. This information is of 
fundamental importance and can be hardly obtained by 
any other means. 

In this letter, we present a calculational scheme that 
allows for easy evaluation of FCS in a multi-terminal 
mesoscopic systems. It is of a great intellectual enjoy- 
ment that this scheme is a simple and universal one. In 
fact, it is hardly more complicated than a conventional 
circuit theory of electric transport and is based on slight 
extension of Kirchoff rules to 2 x 2 matrix structures. 

We start by introducing current operators 7j , each be- 
ing associated with the current to a certain terminal i. 
Extending the method of (hJ we introduce a Keldysh- 
type Green function defined by 

- h + \n £ Xi(t)ii) ® G(t, t') = 5(1 - i') (l) 

i 

Here we follow notations of a comprehensive review [ p"T| , 
Xi are time-dependent parameters, T3 is a 2 x 2 ma- 



trix in Keldysh space, H is the one-particle Hamiltonian 
that incorporates all information about the system lay- 
out, including boundaries, defects and all kinds of elas- 
tic scattering. We use "hat", "bar" and "check" to de- 
note operators in coordinate space, matrices in Keldysh 
space and operators in direct product of these spaces re- 
spectively. The Eq. [l] defines the Green function un- 
ambiguously provided boundary conditions are satisfied: 
G(t,t') = G(x,x';t,t') approaches the common equilib- 
rium Keldysh Green functions G] (t — t') provided x, x' 

are sufficiently far in the terminal i. These g\ (t — t') 
incorporate information about the state of the terminals: 
their voltages Vi and temperatures T^. 

One can easily see by traditional diagrammatic meth- 
ods [JllJ that the expansion of G in Xi(t) generates all 
possible diagrams for higher order correlators of Jj(t) 
and thereby incorporates all the information about statis- 
tics of charge transfer. If we limit our attention to low- 
frequency limit of current correlations, we can keep time- 
independent \i. In this case, the Green functions are 
functions of time difference only and the Eq. 1 separates 
in energy representation. It is convenient to introduce 
the following x^-dependent action defined as a sum of 
closed diagrams: 
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(2) 



This allows us to express the probability for iVj electrons 
to be transferred to the terminal i during time interval 
to 



P({Ni}) 



n 



d Xi -SUxi})-iY..N iX i 



2tt 



(3) 



(Higher-order) derivatives of S with respect to \i gi ye 
(higher-order) moments of P({Ni}). First derivatives 
yield average currents to terminals, second derivatives 
correspond to the noises and noise correlations. 

Using special properties of current operators, x~ 
dependent terms in Eq. [l] can be gauged away. [[To[ |H| 
The x dependence of G is thereby transferred to the 
boundary conditions: the gauged Green function far in 
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FIG. 1: The graph of the circuit theory, associated with a 3 
terminal mesoscopic system. 



each terminal shall approach Gj(e) defined as 

G;(e) = exp(iXif 3 /2)Gf (e) exp(-*x,7s/2) (4) 



In the present form, the Eq. (|lj) with relations (|4|| 
solves the problem of determination of the FCS for any 
arbitrary system layout: one just has to find exact 
quantum-mechanical solution of a Green function prob- 
lem. This is hardly constructive, and we proceed further 
by deriving a simplified semiclassical approach. First, we 
note that even in its exact quantummechanical form the 
Eq. 1 possesses an important property. We consider the 
quantity defined similar to standard definition of current 
density j Q (x,e) = lim x ^ x >(V' a - V Q )G(x, x'; e)/m. By 
virtue of Eq. 1 this quantity conserves so that 



dx a 







(5) 



This looks like the conservation of particles at a given 
energy. However, this relation contains more information 
since it is a conservation law for a 2 x 2 matrix current. 

Next, we construct a theory which makes use of this 
conservation law. We concentrate on the Green function 
in coinciding points, G(x, e) = iG(x, x'; ej/iris. So defined 
Green function has been introduced in several semiclas- 
sical theories. [0 |l^, [l4| It satisfies the normalization 
condition G 2 = 1. We relate the "current density" j 
to gradients and/or changes of G{x), very much like the 
electric current density is related to the voltage in circuit 
theory of electric conductance. Following the approach of 
the circuit theory, we separate a mesoscopic layout into 
elements: nodes and connectors, so that the G(x) is con- 
stant across the nodes and drops across the connectors. 
One may associate a graph with each circuit, so that its 
lines would denote the connectors, and internal and 
external vertices correspond to the nodes and terminals, 
respectively. (See Fig. 1) This separation of actual lay- 
out is rather heuristic, similar to separation of an electric 
conductor of a complicated geometry onto nodes and cir- 
cuit theory elements. The bigger the number and the 
finer the mesh of the nodes and connectors, the better 
the circuit theory approximates the actual layout. The 
nodes are similar to the terminals, the difference is that 
G is fixed in the terminals and yet to be determined in 



the nodes. The G in nodes are determined from Kirchoff 
rules reflecting the conservation law (^): sum of the ma- 
trix currents from the node over all connectors should 
equal zero at each energy. For this, we should be able 
to express the matrix current via each connector as a 
function of two matrices Gij at its ends. 

The connector can be quite generally character- 
ized by a set of transmission eigenvalues T„ .[O, |l5| ] 
The problem to solve is to express matrix current via the 
connector in terms of Gj(j). This problem shall be ad- 
dressed by a more microscopic approach and was solved 
in |l4| for Keldysh-Nambu matrix structure of G. It is 
a good news that the derivation made in (m) does not 
depend on concrete matrix structure and can be used for 
the present problem without any modification yielding 



Z~Y [dE- 



Tk ij) [GuGj] 



+ ({G i ,G j }-2) 



(6) 



Each connector (i, j) in the layout contributes to the total 
Xi-dependent action (|J). The corresponding Sy contri- 
bution reads :fll3| 



Sij(x) 



' ^ 



dETi In 



jT^ ({G, 



Gj} 



(7) 

Now we are ready to present a set of circuit theory 
rules that enables us to evaluate the FCS for an arbi- 
trary mesoscopic layout, i. The layout is separated onto 
terminals, nodes, and connectors, ii. The Gj in each 
terminal j is fixed by relation (^) thus incorporating in- 
formation about voltage, temperature and counting field 
X in each node. iii. For each node k 7 the matrix cur- 
rent conservation yields a Kirchoff equation 1^ = 0, 
where the summation is going over all connectors (i, k) 
attached to node k, and Iik are expressed with (|6|) in 
terms of Gj(M. iv. The solution of resulting equations 
with condition G\ = 1 fixes Gk in each node. v. The 
total action S(x) is obtained by summing up the con- 
tributions Sij({xi}) of individual connectors, those are 
given by (@): S({ Xl }) = &,•({»}) vi. The statis- 

tics of electron transfer is obtained from the relation (0). 



It is time to discuss the limits of applicability of the 
whole scheme. By virtue of semiclassical approach, the 
mesoscopic fluctuations coming from interference of elec- 
trons penetrating different connectors are disregarded. 
So that, we assume that conductivities of all connectors 
are much bigger than conductance quantum e 2 /irh. The 
same condition provides the absence of Coulomb block- 
ade effects in the system. Besides, we have disregarded 
the possible processes of inelastic relaxation in the sys- 
tem. The latter can be eventually taken into account but 
it would considerably complicate the scheme. The point 
is that the inelastic scattering would mix up the G at 
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FIG. 2: The logarithm of the current probabilities in the 3- 
terminal chaotic quantum dot as a function of I3, under con- 
dition I\ = Ii. The insert presents the system configuration. 
The resistances R of all connectors are assumed to be equal. 
1 - tunnel connectors, 2 - diffusive connectors, 3 - ballistic 
connectors. 



eigenvalues in each connector k: Zk(x) = Y2 /[2 + 

n 

can be evaluated for any particular dis- 
tribution p(T) of transmission eigenvalues in the given 
connector and completely defines its scattering proper- 
ties. For a example, if we denote Rq = irh/e 2 , then 
R^ 1 = 2Rq Zf.(l) is an inverse resistance of the con- 
nector. One can also express the Fano factor F = 
(T(l - T))/(T) as F = 1 - 2(d/dx)]ogZ(x)\ x=v With 

the use of Zk(x) the conservation law can be efficiently 

3 

rewritten as [ P k G k , G 4 ] = 0, where p k = Z k (g k ■ g 4 ). 
fc=i 

The latter enables one to look for the vector g 4 in the 
3 

form g 4 = M _1 ^2 p k gk, with M(x) being an unknown 

k=l 

normalization constant. Using the condition g 4 • g 4 = 1 
we obtain the set of equations 

3 3 

v l = Zi (m- 1 9ij (x) /'')• a 2 X 9a (x)p¥ 

(8) 

where 



different energies, so that one can not solve the circuit 
theory equations separately at each energy. 

As an illustration of the presented scheme, we will con- 
sider in the rest of the paper the FCS of the 3-terminal 
chaotic quantum dot. The system is sketched in the in- 
set of Fig. 2. The heuristic circuit, associated with this 
mesoscopic system is shown by dashed lines. It includes 
only 3 external vertices, corresponding to terminals, 3 
arbitrary connectors, associated with the contacts, and 
the node {4}, representing the quantum dot itself. This 
separation is valid provided the cavity is in the quantum 
chaotic regime. (See [ JL6| for definition). This regime 
corresponds to full isotropization of the Green function 
G(x, x' , e) within the dot, so that G 4 (e) can be regarded 
as a constant at a given energy. 

Since the normalization G\ = 1 holds for each ver- 
tex, we use the parametrization Gk — gk ■ t, gk ■ gk = 
1. Here gfc is a 3-D vector, and r = (Ti,f 2 ,f 3 ). In 
the absence of counting fields the Green functions in 
the terminals are given by a zero condition G k = 

I 0/1 \ \ nt^ k -t I , where Fermi distribution func- 

V _ Jk) l Jh - 1 / 

tion fk{E) = {cxp[(£ - eVk)/T k ] + l}" 1 accounts for the 
bias voltages Vk and the temperatures Tk in the termi- 
nals. The Xi-dependence of Gk(x) is then given by Eq. ^. 

We see that G 4 (x) = g 4 • r is in fact the only func- 
tion to be found. For that, we proceed by applying 
the current conservation law, Y], Ik,i = 0, inside the 
dot. We present the currents Ik,A given by (Q) in the 
form IkA = \Zk{gk ■ g 4 )[Gfc,G 4 ], the scalar function 
Zk(x) incorporating the information about transmission 



9H(X) = ft(x)-ft(x) = (1-2/0(1 -Vj) 

+ 2 e *(x*-Xi )/,(!- f 3 ) + 2 e -fe-xi )f j (l-f i ) 

The fixed point {p\(x\ M*(x)} of the mapping (^) finally 
yields G 4 in question. 

The total action can be found by applying the rule (v) 
at page || and reads 

3 3 
S M = / d ^,{M-\x)Y.9Mp l M) (9) 

n t=l J ] = 1 

The partial contributions Si(x) in the above equation 
should be determined from the relation J^Si(x) — 
-Z i {x),S i {l) = Q. 

We specifically consider three particular types of con- 
nectors: tunnel(T), ballistic(B) and diffusive(D). The 
corresponding contributions to action are: St(x) = 
-±{Ro/R)(x - 1), S B (x) = -{Rq/R) log[(l + x)/2], 
S D (x) = -\(R /R)log 2 (x + Vx^T) (Tf|, R being 
a resistance of the connector. For tunnel connector 
T n -c 1 for all n. For ballistic connector N chan- 
nels are opened (T n w 1 for n < N), and the other 
are closed. In the diffusive connector the transmission 
eigenvalues are distributed according to universal law 
p(T) = i? /2i?TVT — f. 

Analytical results for FCS (0) are plausible only for 
the system with tunnel connectors. To assess general sit- 
uation we solve the Eqs. (^) for given Xi numerically. To 
find the probability distribution, we evaluated the inte- 
gral (|3|) in the saddle point approximation, assuming \% 
to be complex numbers. Saddle point approximation is 
always valid in the low frequency limit we consider, since 
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in this case both action S and number of transmitted par- 
ticles Ni — Iito/e ^> 1. Due to the current conservation 
law only two of three counting fields Xi are independent, 
and one can set x 3 = 0. The relevant saddle point of 
the function fl(x) — S(x) + ixihh/e + «X2^2^o/e al- 
ways corresponds to purely imaginary numbers {xt, X2} 
in the upper half plane. The probability reads P{I\, I 2 ) ~ 
exp[— f2(x*)]. Evidently, fi(x*) is the Legendre transform 
of the action, and it can be regarded as implicit function 
on/( X *). 

In the following we assume the short noise regime 
eV ^> kT when the thermal fluctuations can be disre- 
garded. The energy integration in (^) becomes trivial, 
since fi(e) = or 1, and it is sufficient to consider only 
the case Vx = V2 = 0, V3 — V. Any other possible setup 
can be reduced to the number of previous ones by appro- 
priate subdividing a relevant energy strip. The results of 
these calculations are shown in Fig. 2 and 3. We see that 
the maximum of probability occurs at I\ = I2 = — V/3R, 
I3 = 2V/3R that simply reflects the usual Kirchoff rules. 
The current distribution P(Ii,l2) for a ballistic system 
is bounded. It is due to the fact, that Zb{x) contains 
the finite number of open channels, contrary to the tun- 
nel or diffusive type configurations, where it is not the 
case. From (||) and (^) we can also find a zero noise 
and noise correlations matrix $y = eR~ 1 VFij, F%x = 
F 22 = (4 + 3F)/27, F 33 = (4 + 6F)/27, F 12 = -2/27, 
F13 = F 23 = -(2 + 3F)/27, where F is a Fano factor. 
Since Fb = 0, Fo = 1/3 and Ft = 1 one concludes, that 
for a fixed average currents through connectors the Gaus- 
sian's currents fluctuations will increase in the sequence 
ballistic^diffusive— >tunnel. Fig. 2 and 3 show, that the 
similar behavior is also traced in the regime of the large 
current fluctuations. The essential point here is that the 
cross-correlations always persist regardless the concrete 
construction of the connectors. For the case of multilead 
chaotic cavities the results for shot noise in our theory 
coincides with those, obtained by means of random ma- 
trix theory J lj, and with use of "minimal correlation 
approach" 

In conclusion, we present a constructive theory for 
counting statistics for electron transfer in mesoscopic sys- 
tems. With this theory, one can easily make theoretical 
predictions for all FCS, thereby facilitating experimental 
activities in this direction. Up to now, only the noise 
has been measured. In our opinion, the measurements 
of FCS can be easily performed with threshold detectors 
that produce a signal provided the current in a certain 



terminal exceeds the threshold value. If the threshold 
value exceeds the average current, the detector will be 
switched by relatively improbable fluctuations of the cur- 
rent described by FCS, its signal being proportional to 
the probability P{I\, I 2 ) of these fluctuations. 
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FIG. 3: The contour maps of the current distribution log[P(/i, I2)] in the 3-terminal chaotic quantum dot for different config- 
urations of connectors, (a) - ballistic connectors, (b) - diffusive connectors, (c) - tunnel connectors. 



